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ABSTRACT 



Context. In the course of a mid-infrared imaging campaign of close-by active galaxies, we discovered the mid-infrared counterparts 
of bright compact radio sources in the central star-forming region of NGC 1 808. 

Aims. We aim at confirming that these sources are deeply embedded, young star clusters and at deriving some of their intrinsic 
properties. 

Methods. To complement the mid-infrared data, we have collected a set of near-infrared data with ISAAC at the VLT: J, Ks, and L' 
images, as well as low-resolution, long-slit spectra for three of the sources. 

Results. Surprisingly, the new images unveil a near-infrared counterpart for only one of the mid-infrared/radio sources, namely M8 in 
the L' band. All the other sources are so deeply embedded that their emission does not pop out above an extended diffuse near-infrared 
emission. The near-infrared spectra of the sources look alike, with intense, ionised hydrogen lines. This supports the interpretation of 
these sources in terms of embedded young clusters. We derive extinctions and ionising photon production rates for two of the clusters. 

Key words. ISM: dust, extinction, ISM: HII regions, Galaxies: star clusters, Galaxies: individual: NGC1808, Infrared: galaxies 



1. Introduction 

Star clusters appear to be the most frequent mode of star for- 
mation in galaxies. Therefore, the understanding of star cluster 
formation and evolution is necessary for studying the stellar pop- 
ulations in galaxies in general. Young clusters in starburst envi- 
ronments have been extensively studied since the beginning of 
the nineties, thanks to the high angular resolution of the Hubble 
space telescope (HST). Yet, the very early stages of their evo- 
lution are still poorly known, as they are expected to form in 
dense dusty environments, hence to suffer high extinction. As it 
is only recently that high angular resolution instruments in the 
near- and mid-infrared (NIR, MIR) have become available, only 
a small sample of embedded young clusters is known. 

In a previous study dGalliano et al.l I2005I) we discovered 
bright MIR sources within the central starburst of NGC 1808, 
associated with already known radio sources. We suggested that 
these sources were young embedded star clusters. The present 
research note aims at complementing the MIR dataset with new 
NIR data for these sources. We present NIR images of the cen- 
tral starburst, as well as NIR long slit spectra of three of the 
detected MIR/radio sources. The spectra show intense nebular 
emission lines and therefore bring strong support to an interpre- 
tation of the sources in terms of young embedded clusters. From 
the spectral measurements, some basic parameters such as ex- 
tinction and ionising photon emission rate are derived for the 
embedded clusters. 



NGC 1808 i s loc ated at a distance of 10.9 Mpc 
dKoribalski et al.l [l996), which corresponds to a scale of 
53 pc per arcsec. Its redshift is 0.003319 (TKoribalski et al.l 
2004). The central region (inner 750 pc) of this galaxy is 
und ergoing an intense episode of star formation, discussed in 
e.g. lTacconi-Garman et all dl996ll2005l) . 



2. Observations 

With ISAAC at VLT/UT1, we collected NIR images and spec- 
tra of the circumnuc lear MIR sources discovered in NGC 1808 
dGalliano et al.ll2005l). We adopt hereafter the source nomencla- 
ture as in iGalliano et al.l (120051) . We recorded three images: J 
(1 .2 /imfl Ks (2.2 /im) and L' (3.8 /imfl The image angular res- 
olution is around 0.6". We also obtained low resolution long 
slit spectra for two slit positions. For one slit position (pass- 
ing through sources M6 and M8), we obtained a 2 fim spectrum 
(1.8-2.4/mi) and a A fim spectrum (3.1-4.1 /mi). For the other 
slit position (across the nucleus Ml and source M3), we could 
only obtain a 2 fim spectrum. The spectral resolutions are R=450 
and 360 for the 2 fim and 4 fim spectra respectively. The slits 
were positioned performing offsets refe renced on the nucleu s 
Ml, according to the measurements in IGalliano et al.l {2005). 
The data reduction was performed using the ECLIPSE dedicated 
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Fig. 1. The central 18 "xl8 "region of NGC1808 (1 "=53 pc). The axis units are arcsec. Th e nucleus (source Ml ) is positioned at 
coordinates (0"; 0"). North is up and East is to the left. The MIR sources are labelled as in iGalliano et all d2005h on the TIMMI2 
map. Zooms towards the two regions delimited by the squares (Regions A and B) are displayed in Fig. [2] The long rectangles feature 
the spectroscopic slits: small ticks are spaced by 1 ", while long thick ticks highlight the positions at which line emission has been 
recorded. 
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Fig. 2. Zooms towards regions A and B of the NGC1808 central starburst. These regions are defined in Fig.Q] The axis units are 
arcsec, North is up and East is to the left. Regarding the slit drawings, see the caption of Fig.Q] 



routines and IRAF. The spectra were extracted through slit win- 
dows of 1 .4 " along the slit. In addition to these data, we use in 
this research note an archive, pipe-line reduced, WFPC2 nar- 
row filter F658N (655 nm-661 nm) image of the starburst, as 
w ell as the TIMMI2 1 0.4, 11.9 and 12.9 /mi images presented 
in lGalliano etail d2005l) . 



3. Results 

3.1. Images 

Figure Q] shows the central starburst in NGC1808 (18"xl8", 
corresponding to 1 kpcxl kpc) at six different wavelengths from 
visible to radio: 0.66 //m, 1.2 /mi, 2.2/mi, 3.8/mi, 12.9 /mi and 
3.6 cm. The bulk of the emission consists of an 20"xl0" emit- 
ting region roughly aligned along PA=135°. In the maps shown 
in Fig.Q] the coordinates (0"; 0") correspond to the brightest 
source Ml . The nature of this source has not yet been clarified in 
the literature. It is suggested to be a mixture of star formation and 
a low luminosity or fading active g alactic nucleus ( Krabb e~et alj 
1200 U IjTmenez-Bailon et al J 12005!). This source is used for the 
relative registration of the different maps, assuming that the lu- 
minosity peak coincides at all studied wavelengths. The follow- 



ing images are displayed in Fig.Q] in order of increasing wave- 
length: 

- the narrow filter WFPC2 images (0.66 //m), which, at the 
redshift of NGC1808 traces the Ha emission, and is hence 
a good tracer of un-embedded star formation in the region, 

- the J (1.2 /mi) and Ks (2.2 /mi) images which trace the pho- 
tospheric emission of stars, less subject to extinction than in 
the visible, 

- the L' band image (3.8 /mi) which can highlight hot dust, and 
hence globally traces star forming activity and AGN activity, 

- the 12.9 /mi image which traces warm dust emission and 
[Nell] line emission. As for the L' band emission, the im- 
age pinpoints star forming regions and the AGN, 

- the 3.6cm image which corresponds to the free-free emis- 
sion from HII regions, to the diffuse synchrotron emission 
produced by electrons escaping SN remnants and being ac- 
celeration by galactic magnetic fields, and to the likely emis- 
sion from recent SN explosions. 

The maps on Fig.Q]bring new details about the extended star 
forming region. Through its Ha emission (WFPC2 image), the 
region looks extremely patchy, likely the result of a complex dis- 
tribution of star formation "nests" and of molecular/dusty clouds 
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Fig. 3. ISAAC spectra of the MIR/radio sources Ml, M3, M6 and M8 in the central region of NGC1808. The slit positions, as well 
as the source positions along the slits are shown in Fig.Q]and[2] 
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Table 1. Measurements for the clusters Ml, M3, M6 and M8 in NCG1808 

Ml M3 M6 M8 



5 



positions with respect to Ml [arcsec] 

Ao- [arcsec] -1.99+0.1 5.82±0.1 2.71+0.1 

AS [arcsec] -3.76+0.1 -5.25+0.1 -5.51+0.1 



flux densities [mJy] 
value value value value 
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: 4. 
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radio data (Collison et al. 1994) 

or!>™ m -0.61 -0.28 - -0.57 

3:6 cm [mJy] 8.96 0.93 0.58 0.76 

* These flux values and their associated uncertainties must be taken with caution since the atmospheric features are strong in this region of the 
spectrum. 



producing a heavy extinction. It is interesting to note the spi- 
ral arm-like structure of the F658N emission, particularly to the 
South of the nucleus. The J and Ks images do not add much 
information to this picture, except for an enhanced contrast for 
some sources which look faint inn the R image. They do not 
unveil any new source, with respect to the R map. 

In the L' band, the region exhibits a bright diffuse emission, 
possibly the signature of a large amount of dust heated by the 
various star forming nests. Several compact sources appear on 
this map. 

The following two images (12.9 /im and 3.6 cm) show deeply 
embedded sources which do not appear on either the R, J, Ks 
or L' maps. This remark is important, showing that embedded 
star forming activity in s tarbursts can be hidde n even at 4 /im. 
As already pointed out in Gallia ncTet alj d2005l) . the correlation 
between the 12.9 /im map and the 3.6 cm map is excellent. 



Figure [2] offers zooms on two regions of interest (regions A 
and B) both containing MIR/radio sources. In region A, M8 is 
detected in L' only, while M6 is not even detected in L' . None 
of these sources show up in the NIR. Regarding region B, there 
is no detection either in Ks or L' . Strikingly, on the Ks map, the 
location of M3 even corresponds to a dark area in the diffuse 
emission. 

Conversely, we note t hat the K knots detected by 
iTacconi-Garman et al.l d 1996b . and which are interpreted as be- 
ing embedded young stellar clusters too, are not detected on the 
12.9 /im map. 

Hence, there are two types of embedded clusters in 
NGC1808: one bright in the MIR, bright in the radio and faint 
in the NIR, and the other bright in the NIR and faint in the MIR 
and the radio. This dichotomy can simply be explained if we 
consider that MIR/radio clusters are at an earlier evolution stage 
than the sole K clusters: at a few million years of age, the clusters 
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Fig. 4. Position-velocity diagram for Bry along the slit which 
contains the two embedded clusters M6 and M8. The X-axis 
shows the wavelength expressed in relative radial velocity in 
kms" 1 and the Y-axis gives the relative angular position in arc- 
sec. 



are still deeply embedded in their formation material, which is a 
mixture of ionised gas (source of the radio emission and nebular 
emission lines) and dust (source of the MIR continuum). At this 
stage, the stars are so heavily embedded that they do not appear 
even on the K maps. This locally embedding material is even- 
tually expelled, but the clusters still reside within the kpc scale 
dusty star formation region, hence suffer some extinction. At this 
stage, the clusters appear as K extincted sources, with no MIR 
or radio counter parts. This in terpretation is further supported by 
iTacconi-Garman et alJd2005l) : not only is the PAH-to-continuum 
ratio at the K knot positions low, but al so it is high at the ra- 
dio/M IR knot positions (See Fig. 4b in ITacconi-Garman et alj 
120051) . This could mean that much of the PAH/dust around the K 
knots is photoionised, hence weakening the 3.3 yum PAH feature, 
has been photodissociated or has simply been blown away, while 
this has not yet happened in the MIR/radio knots. 

3.2. Spectra 

We have obtained NIR spectra of three embedded clusters and 
the nucleus: 2/im and 4/im spectra for M6 and M8, and only 
one 2 /mi spectrum for Ml and M3. The projected slit positions 
are drawn in the images in Fig.Q]and Fig. [2] On these drawings, 
small ticks mark each arcsec, while the long and thick ticks mark 
the positions at which line emission has been recorded. We de- 
tect emission lines for Ml (the nucleus), M3, M6 and M8, and 
we even detect line emission in the region between M6 and M8 
devoid of any conspicuous MIR or radio source. The spectra ob- 
tained for Ml, M3, M6 and M8 are shown in Fig. [3] 



For Ml, the continuum shown on Fig. [3] is, at the achieved 
angular resolution, associated to the emission line source. Three 
emission lines are un-ambiguously detected and measurable: 
Pace, H2 1-0S(1) and Bry. The Hel 2.058 jum is only tentatively 
detected and its flux cannot be measured because of the pres- 
ence of artificial features on this part of the spectrum. Deep 
2.3 /urn CO absorption bands are detected, tracing the presence of 
red super giants. The equivalent width of Bry is 4.5A. This ca n 
be compared to Starburst99 predictions (Leithe rer et al.l [l999). 
but we must keep in mind that part of the continuum, even in 
K, may arise from hot dust associated to the AGN. Hence, the 
measured value of the equivalent widt h of Bry is a lower limit. 
By comparing with Fig. 89 and 90 of Leith erer et alJ dl999h . we 
deduce that (1) in the instantaneous star formation case, such an 
equivalent width implies that the cluster is younger than 7 Myr; 
(2) in the continuous star formation case, there is no age con- 
straint, since we ignore the AGN 2 fim continuum share. In con- 
clusion, the K Ml spectrum does not clarify the nature of the 
source. 

In the following, we focus on the three circumnuclear 
sources M3, M6 and M8. Their spectra look quite similar in 
shape and are typical HII region spectra. Emission lines from 
atomic (Pa<?, BrS and Bry) and molecular hydrogen show up. 
In M6 and M8, for which we have an additional L' spectrum, 
we detect the hydrogen lines Pay and Bra as well as the PAH 
band at 3.3 fim. This strongly supports an interpretation of the 
MIR/radio source s in terms of heavily e mbedded clusters as al- 
ready suggested in Galli ano et al.l d2005l) . In addition, such spec- 
tra fix an upper limit of 10 Myrs to th e cluster ages, otherw ise no 
nebular gas emission would be seen (lLeithereretalJl!999l) . The 
continuum shown on the spectra of M3, M6 and M8 is not asso- 
ciated with the sources of emission lines. This continuum only 
shows the diffuse emission of the starbursting region. While the 
continuum emission is extended, the emission lines are only ob- 
served at the location of the MIR sources along the slit. Hence, 
the CO absorption bands seen on the spectrum continua are not 
associated to the embedded clusters. 

Fig. [4] displays a 2-D spectrum for the Bry line along the 
M6-M8 direction. The line emission extends over several arcsec 
(~ 250 pc) and exhibits a complex velocity structure, sugges- 
tive of an outflow of ionised material. This probably witnesses 
the expelling of the locally embedding material suggested in 
Sect.O 

3.3. Measurements 

Table Q] provides a summary of the measureme nts performed 
on the data and gives the radio measurements of ICollison et al.l 
d!994h . The following procedure has been followed for the flux 
density measurements and the uncertainty estimates. The rela- 
tive positions of the three sources with respect to Ml are taken 
as those in the TIMMI2 12.9 fim image. They are given in the 
table. At these positions, we measure the flux densities inside 
a 0.6" radius aperture. The background is estimated by com- 
puting the median within an annulus with radii 0.6 " and 1 .2 ". 
Since, for these data, the background is difficult to interpret, we 
consider two values of the flux density for each source, one with 
and one without background subtraction, this is what we call 
the measurement error. We then compute an error defined as the 
quadratic sum of the measurement error and the photometric cal- 
ibration uncertainty. The following flux calibration uncertainties 
are considered: 10% for the R, J and Ks images and 20% for 
the L' and MIR images. The two values presented in table Q] 
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correspond to the final lower and upper limits of each measure- 
ment. When no clear identification of the source in the consid- 
ered aperture is possible (for example, for the data in the R, J 
and Ks band), we did not consider the lower limit defined above, 
but only the upper limit. For the line fluxes, we consider a 20% 
uncertainty, dominated by the calibration. Nevertheless, for the 
lines on the blue side of the 2 /jm spectrum (Paa, Br<5 and H2 1 - 
0S(3)), the uncertainty may be greater than this value because of 
the strong atmospheric features present in this wavelength range 
and the values given must be taken with caution. 

4. Cluster derived parameters 

In this section, we discuss some intrinsic parameters of the 
sources which can be derived from the NIR emission lines. The 
procedure has been applied on M6 and M8, for which both the 
2 yum spectra (Bry flux) and 4 fim spectra (Bra flux) are available. 
Under the assumption of simple foreground extinction, we de- 
rive the extinction towards the line emission gas, and then com- 
pute the un-extincted line luminosities. This in turn allows an 
estimation of the ionising photon emission rate. 

We have deduced the extinction from the Bry/Bror ratio. 
This ratio has a low dependency on temperature and density. 
Assuming reasonable values for the temperature and density 
offfll regions, T e = 10 4 K and N e = 10 4 cm" 3 dOsterbrockl 
119891) . we obtain figures of Av=10.2 and 15.0, for M6 and 
M8 respectively. This leads to un-extincted Bry fluxes of 
8.9xl0~ 15 erg s _1 cirT 2 in M6, and 9.5xl0~ 15 erg s -1 cnT 2 in 
M8. 

For an embedded young star cluster, the ionisation-bounded 
case is very likely: the HII region is optically thick in the Lyman 
continuum and all the stellar ionising photons are absorbed. 
Then, the number of hydrogen ionising photons (A < 912A) is 
directly proportional to the flux in any specific recombination 
line. 

From lOsterbrockl d 19891) . 
Q[H + ] = X £22L , where aa/afi, ~ 2.96. (1) 

In case B, we getL,Ha7LBy= 103.6. The ionising photon emission 
rates of the two embedded clusters can then be derived and are 
found to be around 9.6xl0 51 s _I and 10.2xl0 51 s~\ in M6 and 
M8 respectively. This corresponds to about 1000 06 stars in each 
cluster. 

The way we have treated extinction in the above computa- 
tion is a crude one. This is undoubtly a first order estimate which 
calls for a future, more realistic approach, once additional data 
constraints will become available. Yet, the interpretation of the 
MIR/radio sources in terms of embedded young clusters is con- 
firmed and, in addition, an upper limit on the cluster ages is ob- 
tained at lOMyrs. 



- The fact that the MIR/radio sources do not show up in the 
NIR images demonstrates that such embedded clusters must 
be searched for primarily in the MIR or at radio wavelengths. 

- A complex velocity structure in the ionised gas is observed 
around sources M6 and M8, suggesting the presence of clus- 
ter winds. 

- For two sources, M6 and M8, the extinction towards the neb- 
ular gas has been estimated, with values of Av=10.2 in M6 
andAv=15in M8. 

- The extinction-corrected Bry fluxes for M6 and M8 imply 
ionising photon emission rates of respectively 9.6 xlO 51 s _1 
and 10.2xl0 51 s _1 . 

A more complete analysis of the deeply embedded clus- 
ter parameters, in particular the derivation of their age and 
mass, requires some knowledge of their MIR spectral features 
dGalliano et al.|[2008h and it is left for a subsequent analysis, af- 
ter spectral MIR data will have been collected. 
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5. Conclusions 

The main conclusions are: 

- None of the MIR/radio sources unveiled by iGalliano et alj 
(2005) can be detected in the J, Ks and even L' band images, 
except for M8 in the L' band. 

- For three of the sources, namely M3, M6 and M8, we have 
collected NIR spectra: they display intense nebular emission 
lines, confirming that the sources are young embedded clus- 
ters. 



